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Abstract 
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namics for an oppositely charged pair of particles, with a view to extending these 
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I. Solving Maxwell's Equations with Dual Potentials 

Dirac^ showed that Maxwell's equations could be extended to include both elec- 
trically and magnetically charged particles by connecting the magnetically charged 
particles to strings. In the absence of magnetically charged particles one can ap- 
ply Dirac's method to ordinary electrodynamics by connecting electrically charged 
particles to strings. In this formulation Maxwell's equations become equations for 
dual potentials C M whose sources are the polarization currents produced by the Dirac 
strings. The potentials themselves depend upon the location of the strings but they 
yield the same string independent electromagnetic fields as the usual procedure. 

If in addition the dual potentials are minimally coupled to Higgs fields, these 
fields necessarily carry magnetic charge. Such a theory describes the motion of elec- 
trically charged particles connected by Dirac strings in a dual superconductor. If 
extended to non-Abelian gauge theory' 2 ! it becomes a concrete realization of the 
Mandelstam 'tHooft' 3 ' 4 ! picture of color confinement as a manifestation of dual su- 
perconductivity. We have found that in order to understand this mechanism for 
confinement it is very helpful to first have a clear picture of how dual potentials 
work in ordinary electrodynamics. Therefore, in this paper we present an elementary 
discussion of Dirac's method applied to electrodynamics and work out some simple 
examples. Our discussion is entirely pedagogical and contains nothing new although 
some specific results obtained here may not be readily accessible elsewhere. Our goal 
here is to make the following paper as comprehensible as ! possible. 

Consider first a sourceless linear dielectric medium. Then Maxwell's equations 

(a) (b) (c) 

V-/3 = V-5 = D = eE (1.1) 

V x H-d D = Vxi + 8qB = B = fiH , 

can be solved by introducing vector potentials in either of two ways. The conventional 
choice is to write 

BeVxA, E = -d A-VA , (1.2a) 
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in which case eqs. (1.1b) become kinematical identities and the dynamics is contained 
in eqs. (1.1a). The vector = (Aq, A) is called the vector potential. The alternate 
(dual) choice is to write 

D = -V x C, H = -d C - VC , (1.26) 

in which case eqs. (1.1a) are kinematical identities and eqs. (1.1b) contain the 
dynamics. The vector = (Co, C) is called the dual vector potential. 

Let us first use to solve the source free Maxwell eqs. (1.1) in order to get 
accustomed to using the dual potential. We first write eqs. (1.2b) in covariant form 
by defining: 

Gok = H k , dj = eijkDk, (1-3) 
so that eqs. (1.2b) take the form 

G lu , = d ll C v -d v C ll . (1.4) 

Then using eq. (1.3) we can write the constitutive 

^ e ijkGjk, Bi = fiGoi, (1.5) 

PljGap = 0. (1.6) 

Eqs. (1.6) for have the same form as the usual Maxwell equations for A^, obtained 
from (1.1a), with the replacement /i — > e, and they are solved in the same way. Eq. 
(1.5) then gives the electromagnetic fields E and B in terms of C M . 

Electric current sources = (p, j) appear only in eqs. (1.1a) and not in eqs. 
(1.1b). Hence in the presence of electric currents eqs. (1.1b) remain valid and are 
still kinematic identities in terms of A^. Eqs. (1.1a), in contrast, are no longer 



In a relativistic medium e = -k 

A* 

equations as 

Ei = 

and Maxwell's eqs. (1.1b) as 
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identities in terms of C M . However, Dirac has shown how to generalize eqs. (1.2b) in 
order to satisfy eqs. (1.1a) with dual potentials even in the presence of electric 
currents. 

When charged particles are present eqs. (1.1a) become 
Gauss' Law Ampere's Law 

- - - «n (1-7) 

V-D = p, Vxff = J + f. 

Suppose that the total charge Q = J pdx = 0. (If Q ^ 0, then there will be Dirac 
strings extending to infinity, but nothing essential will be changed.) Then we can 

— * — * 

always find a polarization vector P and a magnetization vector M so that 

BP 

p = -V-P, j=VxM+— , (1.8) 

— * — * 

where P is the dipole moment per unit volume and M is the magnetic moment per 
unit volume. Inserting eq. (1.8) into (1.7), we obtain 

V • (D + P) = 0, V x (H - M) - ^ + P) = 0. 

ot 



Hence, 



dC 

D = -V xC-P, H = -VC - — + M. (1.9) 

ot 



Eqs. (1.7) then become kinematical identities and eqs. (1.1b) contain the dynamics 
as before. Using the definitions (1.3) of G^, we can write eqs. (1.9) in the covariant 
form 

G^ = d^C v -d v C^ + G s ^ (1.10) 
where the tensor G s ^ v has components 

G s ok = M k , Gfj = -e ijk P k , (1.11) 

and we have now specialized to the case where P and M arise from Dirac strings 
connecting the charged particles; hence the superscript "s" on G s ^ v . Eq. (1.10) is just 
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the generalization of eq. (1.4) to account for the presence of charged particles. Eq. 
(1.11) shows that G s ^ v is the dual of the polarization tensor. Eqs. (1.1b) and (1.1c) 
are unchanged so that eqs. (1.5) and (1.6) remain the same as does the definition (1.3) 
of G^y. The effect of the charged particles is to change the relation (1.4) between G^ v 
and to (1.10) where is determined in terms of p and j by solving eqs. (1.8) 

— * — * 

for P and M. This is done in Section II for a pair of oppositely charged particles. 
Substituting eq. (1.10) into eq. (1.6) we obtain 

d a p(d a C p - d p C a ) = -d a piG s ap , (1.12) 

which determines the dual potentials in terms of G s ^ v . Eqs. (1.12) provide an 
alternate form of Maxwell's equations which are completely equivalent to the usual 
form expressed in terms of the vector potential A^, namely 

d a e{d a Ap-d p A a ) =jp. (1.13) 

All text books on electricity and magnetism could be rewritten using only dual po- 
tentials C M satisfying eq. (1.12) and the same electromagnetic forces between charged 
particles would be obtained. The potentials themselves, however, could be completely 
different. For example in a dielectric medium having a wave number dependent dielec- 
tric constant e(q) — > as q 2 — > (corresponding to antiscreening at large distances), 
the potentials A^ determined from eq. (1.13) would be singular at large distances, 
while the dual potentials satisfying eq. (1.12) with /i=l/e^ooasg 2 ^0 would 
be screened at large distances. Use of the potentials A^ to describe this system would 
introduce singularities which do not appear in the dual potentials C^. Hence the dual 
potentials are the natural choice to describe a medium with long range antiscreening. 

Note that for /x = e= 1,B — H, D — E and substituting eqs. (1.9) in (1.1b) 
gives the equation for the dual potentials in three-dimensional notation: 

V • (— VCq — 3qC) — —V • M, (1.14) 
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V x (-V x C) + <9 (-<9 C - VC ) = V x P - —. (1.15) 
These equations are identical to eq. (1.12) with fj, — 1. They have the same form as 

— * 

the equations for A^, the ordinary vector potentials in a polarizable medium with P 

—* — * — * 

and M interchanged. For example —V • M is the source of Co. However eqs. (1.14) 
and (1.15) describe the electrodynamics of electrically charged particles moving in 

— * — * 

the vacuum and P and M are the polarization and magnetization respectively of the 
Dirac strings attached to these particles, as we shall now see. 

II. The Fields of a Pair of Oppositely Charged Particles 

We now apply the results of the previous section to the case of two particles of 
charge e(— e) moving along trajectories x\(t)(x 2 (t)) in free space with fj, — e — 1. 
Then 

p(x,t) = e[5 3 (x-x 1 (t))-5 3 (x-x 2 (t))}, (2-1) 

and 

Mt)=e[v 1 5 3 (x-x 1 (t))-v 2 5 3 (x-x 2 (t))}, (2.2) 

where V{ = i = 1,2. We must find a polarization P and magnetization M 
satisfying eq. (1.8) with p and j given by eqs. (2.1) and (2.2). The solution of this 
problem was given by Dirac.' 1 ' Let y((J,t) be any line L(t) connecting x 2 {t) and x\(t) 
i.e., y(ai,t) = xi(t), y(a 2 ,t) = x 2 (t),a 2 < a < a\. (See Figure 1.) On each element 
dy of L place a dipole moment dp = edy. It is evident from Fig. 1 that the charge 
and current density produced by the sum of these dipoles is that due to the pair of 
moving oppositely charged particles, namely eqs. (2.1) and (2.2). To obtain (2.1) 
formally we note that the dipole moment per unit volume P is 

Xl(t) ai 

P(x)=e J dy5(x-$=e jda^^6(x-y(a,t)). (2.3) 

X 2 (t) °2 
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Then 



Xl(t) 



-V-P = -e J dy-V x 5(x-y)=p(x). 

x 2 (t) 



(2.4) 



Furthermore since the line element dy is moving with velocity v = J^2/(o", t), the string 
L in Fig. 1 has a magnetization 



M = ej dyx^5(x-y)=e J da^- x ^S(x - y(a,t)). (2.5) 



Jr-2 



<T2 



— * — # 

Next we show explicitly that eq. (2.3) for P and (2.5) for M give via eq. (1.8) 
the current density. From eq. (2.3) we have 



dP Of 

!n =e di J dyi^tM* - vi^t)) 

x 2 (t) 



ij±6(x - xi(t)) - ^-S(x-x 2 (t)) 



+ 



St 



Xl(t) Xl 

J dy{a, t + 5t)S(x - y(a, t + 5t)) - J dy(a, t)5(x - y(a, t)) 



x 2 (t) 



X 2 



j(x,t) + j t f dy5(x-y). 



(2.6) 



The first term on the right-hand side of eq. (2.6) arises from differentiating with 
respect to the end points with the path fixed. The line integral in eq. (2.6) is over a 
closed contour running from X2(t) to x\(t) along the path y(a,t + 8t) and returning 
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to X2(t) along y(a,t). (See Fig. 2.) We denote y(a, t + dt) — y(a,t) = 5y and the 
element of area dy x 5y = dS. Then by Stokes' theorem 

j t j dy5{x-y) = ~ J dS x V y 5(x-y) 



= e J dyx-^x V x 5(x — y) = — V x M. (2.7) 

Eqs. (2.6) and (2.7) yield eq. (1.8) with j given by eq. (2.2) as asserted. 

Eqs. (1.14) and (1.15) with P and M given by eqs. (2.3) and (2.5) respectively 
determine C^. To obtain the explicit form of the covariant version of these equations 
we note that 



Tl U\ 



f f dy a dy 13 
G'iiv = / dr / da — — 5 (x-y), (2.8) 



T2 (72 



where 



and 



x fi = {t,x), y» = (y ,y), 



dr = yjid^f - {dyf. 

Eq. (2.8) is the standard covariant form for the Dirac string field To show 

that the expressions (1.11) and (2.8) for G s ^ v are the same first set \i = and v — k 
in eq. (2.8): 

Tl 0\ 

/[ dy m dy n 

T2 (72 



7(f *!)/<*-«>=*■ 

ff2 
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Next set fj, — i and v — j in eq. (2.8): 



/[ dy dy 
dr / dr da s3 (x-y) s (yo-t) 

n oi 

/uy 
da—5 3 {x - y) = -e ijk P k . (2.9b) 

Thus eq. (2.8) is just the covariant version of eqs. (2.3) and (2.5). 

Eq. (1.12) with fj, — 1, and given by eq. (2.8) is the covariant form of eqs. 

(1.14) and (1.15) determining the dual potential produced by a pair of oppositely 
charged particles moving in the vacuum. The resulting will depend upon the 
location of the string, but this dependence will drop out in the expression for the 
electromagnetic field tensor G^. We will show in the next section how eqs. (1.14) 
and (1.15) with P and M given by eqs. (2.3) and (2.5) produce the usual expressions 
for the electric and magnetic fields of slowly moving particles. 

To conclude this section we note that the equation of motion (1.12) with /j, — 1, 
namely 

d"G> = 0, (2.10) 
can be obtained from a Lagrangian density C given by 

C = -\G^ = 1 -{H 2 -D 2 ). (2.11) 

We will see in Section V that this Lagrangian gives not only the field equations but 
also the particle equations of motion. 

III. Coulomb's Law and the Biot Savart Law 

To understand better how dual potentials work we will solve eqs. (1-14) and 

(1.15) for slowly moving particles. First consider charges at rest. Then M = and P 
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is time independent, and eq. (1.14) becomes V 2 Cq = 0, (i.e., Cq = 0), and eq. (1.15) 
reduces to 

Vx(-VxC J) ) = VxP, (3.1) 

where we have denoted the static solution C = Cjj (for Dirac). Eq. (3.1) has the 
form of the equation for the vector potential due to a polarization current produced 
by the superposition (2.3) of dipoles. Thus, Cjj is just the vector potential produced 
by a superposition of point dipoles of strength — edy distributed uniformly along the 
string (see eq. 2.3), i.e., is given byt 5 ' 



Xl ^ ^ 

d D (a!) = -± [ dyxj^^. (3.2) 
4vr J \x — y\ 6 

X2 



Then 



where 



so that 



Xl Xl 

-VxC„4l [i,-* 1 *-* 



4:71 ) J \x — y| 3 J \x — y\ 3 

X2 X2 



{Xl 
J dy^5{x-y) + ^ 
X2 



X - Xl) [X-X2 



— xi r \x — X2\ 



= P + D C , (3.3) 



6 ( (yX — X\ ) {x — \ 
Dc = Dcoulomb = T~ I - "p ' ( 3 - 4 ) 

47T \ I - Xll' 3 X — X2 



D = — V x Cjj — P = Dc- (3.5) 



The above elementary derivation of Coulomb's law indicates that it really isn't too 
much harder to work with dual potentials and strings than to work with ordinary 
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potentials and localized charges. The string cancellation mechanism in eq. (3.5) is 
depicted in Fig. 3 in which we have taken the string to be a straight line connecting 
X2 and x\. We see that —V x Cjj gives a divergence free field distribution. The 
singular field passing through the line L is cancelled by the singular polarization P, 
leaving a Coulomb field with a source at x\ and a sink at xi- 

Next, let us solve eqs. (1.14) and (1.15) to first order in v\ and V2 and to zero 
order in the accelerations v\ and V2- First look at eq. (1.15). We choose the gauge 
V • C = 0. (Note V • C D = 0). Then eq. (1.14) becomes 

-V 2 C QD = -V • M, (3.6) 

where we have denoted the solution Cq = Cqd- Eq. (3.6) has the form of the equation 
for the scalar potential due to a polarization charge produced by the superposition 
(2.5) of dipoles. Hence Cqd is just the scalar potential produced by a superposition 
of point dipoles of strength edy x y distributed uniformly along the string, i.e., Cqd 

is [6] 

Xl(t) 

Cod = ^- [ (dyx^).^^. (3.7) 
An J \x — y\ 6 

x 2 (t) 

From eqs. (2.5) and (3.7) we see that the time derivative of M and Cqd do not 
contain terms linear in the velocities. The same is true for 8qCd calculated from eq. 
(3.2) with x\ — > xi(t),X2 — > X2(t). Hence to first order in the velocities eq. (1.15) 
reduces to eq. (3.1), and so C = Cjj an d D = Dq. 

— » — * — * 

To calculate H we use eq. (1.9) with C = Co and Cq = Cqd- We first calculate 

Xl(t) 

d e d f (x — y) /N 

x 2 (t) 

The evaluation of the right-hand side of eq. (3.8) parallels that of eq. (2.6) and we 
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obtain 

where the term 

= Hjiiot Savart = |£_* l(t) |3 " X | f _ ^ |3 j > (3-10) 

arises from time differentiation of xi(t) and ^(i) in eq. (3.8) leaving the path fixed. 
The line integral in (3.9), over the same contour occuring in eq. (2.6), arises from 
moving the string keeping the endpoints fixed. Paralleling eq. (2.7) we then apply 
Stokes' theorem to obtain 

e_ f dyx (x-y) _e_ f (dS x V y ) (x - y) 
St J An \x — y\ 3 St J 4ir \x — y\ 3 

Xl{t) 

"/(*"S(*-£#) 

x 2 (t) 

47r J ot \x — y\° 

= -M + VC 0D . (3.11) 

Eqs. (3.9) and (3.11) then yield 

H = - VC 0D + M = H BS . (3.12) 

Thus we see that the Biot Savart magnetic field Hbs comes from the time derivative 
of the limits X2{t) and x\{t) in the integral for Co, eq. (3.8). The remaining string 
dependent part of ^f- cancels the contribution to H coming from Cqjj and M. 
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We now use the solutions for Cjj, Cod, Dq and Hbs to eliminate the fields in the 
Lagrangian L given by 



dx(H 2 -D z ), 



(3.13) 



to second order in the velocities of the charged particles. The Lagrangian L, defined 
as the integral over the Lagrangian density (2.11) then becomes a function L = 
L(x\, X2, v\, V2) only of the positions and velocities of the charged particles. To higher 
order in the velocities one cannot eliminate the field degrees of freedom in L because 
of the presence of radiation. 

For particles at rest we have C = Cp, Cq = 0, H = 0, D = Dq and 



L(v\ = v 2 = 0) 



dx-D% = - — r—^— 
2 6 4?r £1 - 



-> 1 ' 



(3.14) 



where the self-energy has been subtracted. To first order in the velocities, C = Co 

— * — * 

and D = Dq given by eq. (3.4) with x\ — * x\(t),X2 — > X2(t). In other words the 
static field configuration follows adiabatically the motion of the charged particles. 
Furthermore since the Lagrangian L is stationary about static solutions of the field 
equations we have 



/ 



dx-D 2 
2 



4tt|£i(*) -x 2 (t)\ 

valid to second order in the velocities v\ and V2- 

All the velocity dependence in L then comes from J H 2 which to second order in 
the velocities is 



\ J dxH 2 = \ J dx(H BS ) 2 = -\ 



2 AuR 



V1-V2 + 



v\ ■ Rv2 ■ R 



(3.15) 



where the self-energies have again been subtracted out and where R = x\{t) — X2(t). 
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Hence, we obtain for the second order Lagrangian L (first obtained by Darwin t 7 !) 



L{x\,X2\vi,V2) = | / dx[H BS -D c ] 



— e 



AnR 2 47ri? 



v\ • Rv2 • R 

w 



(3.16) 



As a final remark we connect the notation of this paper to that used in our 
previous work! 2 ! on QCD where we have introduced string fields D s and H s defined 

as 



-P, H s = M, 



(3.17) 



so that eqs. (1.9) take the form 



dC 

D = -VxC + D s , H = -VC - — + H S . 

at 



(3.18) 



The fields D s and H s then cancel the string contributions to —V x C and — VCo — ^ 
yielding fields D and H free of string singularities. For slowly moving particles this 
mechanism is explicitly exhibited by eqs. (3.3) and (3.11). 

IV. The Dual Lagrangian and the Equations of Motion of the String' 8 ' 

The action S describing the electromagnetic interactions of a particle of charge e 
and mass mi with a particle of charge — e and mass rri2 is 



S = 



-mi J yj 1 - v\dt - m 2 J \J 1 - v\dt + J d 4 xC, (4.1) 

t 2 t 2 



where C is given by eqs. (2.11) and (1.10). Varying in the action S gives the 
field equation (2.10). To obtain the equations of motion for the particles and for 
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the string which connects them we vary the string coordinates: — > + 5y^ and 

correspondingly vary the particle positions x\ — > a;^ + (fa^ , a;^ — *• ^2 + sucn that 

5^(r) = 5y X (a 1 ,r) 

6x^r)=5y X (a 2 ,r). (4.2) 



Denote 



Then 



where 



dy X , x dy x , 



y > -f - = r- ( 4 -3) 



9o- " ' <9r 



5 (-1 I d^xCTG^j =- l -j d 4 xG^(x)5G s ^(x), (4.4) 



n en 



5G% = -ee^ X(T J dr J da{[5y%' a + y x 5y' a ]5(x-y)+y x y >a d yP 5(x-y)5y' 3 }. (4.5) 
Hence, 



T2 (72 



n (7i 



- | dxG^5G s ^ = e -e, vXa j dr j daG^(y)(5y x y>« + y x 5y' a ) +y x y' a d p G^(yW 



T2 U'2 



Tl (71 , 

6 e^ a J dr J dal - 5y x -^-(G^y fa ) - 5y a ^(G^y x ) 



2 

T2 (72 



+5y f3 y x y' a dpG^ + ^-(5y x G^y' a ) + 4z(8y a G^y M 



dr da 



Tl <7l 



- Jdr J dal -5y X y' a fd p G^ -5y a y X y' P d p G^ 



T2 (72 
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+5yf } y x y' a dpG^ + ^-(5y a G^y x ) 



Tl d\ 



Ijdrj da\ (5y X y' a yP)(-e^ Xa d p G^ -e^ x d a G^ 



T~2 <J2 



d 



+e^ a d x Gn + e ^ Xa j-{5y a G^y x ) }, 



where we have used the fact that the variations of Sy x vanish at T2 and 
Next we use the identity 

-{^tiv\ a d(3 + e^/3A<9 a + e tlvo ,pd\)G lil ' = 2e^\ a pd v G vl1 

and obtain 

Tl 0~1 

- l - J dxG^SG^ = zjdrj da5y x y' a yPe^ ap d v G^ 

T~2 &2 

Tl 

+§ VAa j dT{5xlG^{ Xl )x x - 6x%G^(x 2 )x x ). 

T2 

We must add to the above variation that of the particle action Sp: 



Tl 



SS P = 5 / dr 



-mi \l x\ — mi J x\ 



T 2 



Tl 

J dT(—mixi a 5xi — vri2X2a^x%) 



T 2 
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= J(mix la Sxi + m 2 X2 a 5x2). (4.9) 



T'2 



The total change in the action 5S due to a change in particle coordinates is 

SS = 5S P + 5 (~ I dxG^G^ 



Tl 

dr 



T 2 



5x1 (rnix la + ^e^ Xa G^ {x^x^ 



+5x2 {rn 2 X2a - ^e llv \ ct G^{x2)x^j 



(4.10) 



In proceeding from (4.8) to (4.10) we used the field equation (2.10) which eliminates 
the string contribution to the variation of the action in eq. (4.8). 

If we had introduced further interactions^ of the C M so that d^G^ v ^ 0, then 
there would have been additional variations of the action arising from the first term 
in eq. (4.8). In that case Hamilton's principle 5S = gives 

(71 

e J day' a y' P ^Xapd v G^{ y ) = 0, (4.11) 

in addition to the Lorentz force equations 

mixia = ^e atlvX G^(xi)x\, 



and 



m 2 x 2a = ^e afM/X G^(x 2 )xi (4.12) 



following from eq. (4.10). Thus eq. (4.11) provides a boundary condition along the 
strings upon the current d v G v ^. 
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V. Conclusions 

Our purpose in writing this paper is pedagogical, though of course, motivated 
by our interest in using dual potentials in QCD. We have seen how normal classical 
electrodynamics can be handled completely in terms of dual potentials, and that 
the use of these potentials gives the solution to the conventional Maxwell equations 
for the electric and magnetic fields and leads to the usual Lorentz force law for the 
motion of charged particles. While dual potentials provide a somewhat awkward way 
to solve electrodynamics when charges are present, they are nevertheless the natural 
variables to describe a dielectric medium with long distance ant i- screening. This is 
the underlying reason for their utility in describing long range QCD. 
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Figure Captions 

Fig.l Dirac String L(t) connecting oppositely charged particles. 

Fig. 2 Closed contour describing path of line integral on right-hand side of eq. (2.6). 

Fig. 3 Diagram representing string cancellation mechanism of eq. (3.5). 
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